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Abstract. Calculations of the properties of rotating strange stars with crusts are performed within the framework 
of general relativity. We employ an equation of state (EOS) of strange quark matter based on the MIT Bag Model 
with massive strange quarks and lowest order QCD interactions. The crust is described by the BPS equation of 
state. A significant increase of the stellar radius is found close to the Keplerian (mass- shedding limit) configuration. 
This leads to the disappearance of the gap between the stellar surface and the innermost stable circular orbit 
(ISCO) at very high rotation rates, for a rather broad range of stellar masses. The Keplerian configuration for the 
strange star with crust corresponds to values of J, T/W, Pisco ~ l/^^isco which are about 20% smaller than in 
the case of bare strange stars. Because the Keplerian configuration is achieved due to the increase of the stellar 
oblateness, the Keplerian frequency (of the rotation) remains almost unaltered. The lack of the gap close to the 
Keplerian rotation could imply a more stringent limit on i^isco, if the existence of such a gap is supported by 
observations, as in the source 4U 1820-30 with the upper QPO frequency 1.07 kHz. If such a constraint is taken 
into account (mandatory existence of a gap) the minimum j^isco is about 1 kHz even with the extreme fine tuning 
of strange quark matter parameters. The minimum i^sco is then obtained for the non-rotating configuration with 
maximum allowable mass. The maximum frequency in the stable circular orbit around the strange star with a 
crust is smaller by about 100 Hz than in the case of a bare strange star. During the spin-down of a magnetized 
strange quark star with crust, the crust matter is absorbed in the equatorial region by the strange matter core. 
The deconfinement of absorbed crust matter is a strongly exothermic process, which would infiuence the cooling 
curve of this compact object. 
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1. Introduction 

Since the mid-1980, when these exotic compact stars were 
introduced (Witten 1984; first detailed models in Haensel 
et al. 1986, Alcock et al. 1986), the existence of strange 
quark stars in the Universe remains a matter of lively de- 
bate (for a review, see Madsen 2000). Recently, these hy- 
pothetical objects were invoked in the context of mod- 
eling compact X-ray and gamma-ray sources (Bombaci 
1997; Cheng et al. 1998; Dai & Lu 1998; Li et al. 
1999). Simultaneously, several studies, in which exact 
2-D general-relativistic calculations of stationary rigidly 
rotating configurations were performed, focused on the 
substantial difference in rapid rotation of strange quark 
stars and (normal) neutron stars. This feature stems from 
the basic difference in equations of state of these two 
types of compact superdense stars (Gourgoulhon et al. 
1999; Stergioulas et al. 1999; Zdunik et al. 2000; Gondek- 
Rosinska et al. 2000). In particular, the difference in the 
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distribution of matter in the rapidly rotating models was 
shown to imply significant difference in the outer space- 
time. Specifically, qualitative differences in the properties 
of the innermost stable circular orbit (ISCO) around these 
two classes of rapidly rotating compact objects have been 
found (Stergioulas et al. 1999; Zdunik et al. 2000). The 
comparison of numerical results for the ISCO (in particu- 
lar, the rotation frequency of plasma blobs on the ISCO, 
i^isco) with observations of the kHz quasi-periodic oscil- 
lations (QPO) in some of the low-mass X-ray binaries 
(LMXB) was used to put constraints on the strange quark 
matter (SQM) models. 

A strange star could be covered by a crust of nor- 
mal matter, a possibility which is particularly natural 
in the case of the LMXB (or of a strange star which 
went through a LMXB epoch). The problem of formation 
and structure of a crust on an accreting strange star was 
studied by Haensel & Zdunik (1991) (see also Miralda- 
Escude et al. 1990). Because of its low mass, typically 
^ 10~^ Mq, the effect of the crust on the exterior space- 



2 



J.L. Zdunik et al.: The crust of rotating strange stars 



time is negligible. However, due to its thickness, the crust 
determines the location of the stellar surface. As it has 
been shown in Zdunik et al. (2000), the presence of the 
crust could play an important role in strange stars ro- 
tating close to Keplerian (mass-shedding limit) frequency. 
Namely, the crust undergoes substantial inflation in the 
equatorial plane, leading to the disappearance of the gap 
between the ISCO and stellar surface. The crust mass, 
thickness and other macroscopic parameters depend only 
on the gravitational force exerted by the SQM core (which 
is determined by the core mass and radius), and on rota- 
tion frequency. 

Preliminary, approximate studies of the crust on ro- 
tating strange stars have been performed by Glendenning 
& Weber (1992), in the framework of Hartle's (1967) 
slow rotation approximation supplemented by the self- 
consistency conditions allowing for determination of the 
stellar parameters close to the Keplerian frequencies 
z/rot — vk- The results of this approach were exploited 
by Zdunik et al. (2000). They have shown that at some 
fixed total baryon mass, Mb, of the strange star, the ro- 
tational increase of the equatorial radius, relative to the 
nonrotating case, is quadratic in t'rot/^'K- As Zdunik et al. 
(2000) have found, the quadratic approximation described 
very precisely numerical results of Glendenning & Weber 

(1992), up to I'rot = vk- 

In the present paper we present results of exact (i.e. 
not based on the slow rotation approximation) 2-D calcu- 
lations of the structure of the crust on rotating strange 
star. We show that the quadratic approximation for 
i?oq(t'rot) - i?(0), AfB,crust(t'rot) " MB,crust(0) is vcry pre- 
cise for i/pot ^ 500 Hz, becomes less precise with a further 
increase of Vmt, and badly underestimates the increase of 
i?cq and Afe, crust for v^ot — ^k- Consequently, the pres- 
ence of the solid crust substantially changes the parame- 
ters of the Keplerian configuration, and the properties of 
the ISCO for v^^t — i^k- 

Strange quark stars with a crust could possess a strong 
magnetic field, and their magnetosphere may be similar 
to that of a neutron star. Therefore, their rapid rotation 
might be expected to be slowed down by magnetic braking. 
Slowing down of rotation would imply a compression of the 
crust matter in the equatorial band. As the density at the 
crust bottom cannot exceed the critical value pb at which 
crust matter becomes unstable with respect to irreversible 
absorption by the SQM (due to neutron drip or quan- 
tum tunneling of nuclei through the Coulomb barrier) , the 
slowing down of rotation is accompanied by a decrease in 
the baryon mass of the crust Mb, crust, and a heating due 
to energy release accompanying absorption and deconfine- 
ment of baryons (deconfinement heating). We derive ana- 
lytical formulae expressing the rate of deconfinement heat- 
ing as a function of i>rot and z/rot- Deconfinement heating, 
as well as heating due to non-equilibrium weak interaction 
processes in SQM (Cheng & Dai 1996) will influence the 
cooling curve of rotating magnetized strange quark stars 
with a crust (Yuan & Zhang 1999). 



In Sect. 2 we briefly present the equation of state 
SQM and the method for constructing relativistic 2-D 
models of stationary rotating strange quark stars with 
crust. Numerical results for the parameters of rapidly ro- 
tating strange stars, including the ISCOs, are presented in 
Sect. 3. Exact numerical results for the crust on rotating 
strange stars are fitted using simple, convenient, analyt- 
ical expressions. The deconfinement heating of spinning- 
down strange quark stars with crusts is studied in Sect. 4. 
Finally, Sect. 5 contains discussion of our results and con- 
clusions. 

2. Equation of state and calculation of rotating 
strange star models 

2.1. Equation of state of strange quark matter 

The differences between the density profiles of a strange 
star and a neutron star result from the basic difference in 
the EOS of their interiors. In the case of a rapidly rotating 
compact object (situation relevant to LMXB), the differ- 
ences between the matter distributions within a neutron 
star and a strange star may be expected to imply differ- 
ences in the shape of rotating objects and in the spacetime 
exterior to the star, and in particular, differences in the 
properties of the ISCO. 

Our EOS of strange matter, composed of massless u, 
d quarks, and massive s quarks, is based on the MIT Bag 
Model. It involves three basic parameters: the bag con- 
stant, B, the mass of the strange quarks, rUs, and the 
QCD coupling constant, ac (Farhi & Jaffe 1984, Haensel 
et al. 1986, Alcock et al. 1986). Our basic EOS corre- 
sponds to standard values of the Bag Model parameters for 



strange matter: B — 56 MeV/fm^ 



200 MeV/c^, and 



ac = 0.2 (Farhi & Jaffe 1984, Haensel et al. 1986, Alcock et 
al. 1986). This EOS of strange quark matter will be here- 
after referred to as the SQMl one. It is the same as that 
used i n Zdunik et al. (2000) and Zdunik & Gourgoulhon 
(2001). It yields energy per unit baryon number at zero 
pressure Eq = 918.8 MeV < Ei^'^Fe) = 930.4 MeV. For 
the SQMl EOS the maximum allowable mass for static 
strange stars is M^^^ — 1.8 Mq. 



2.2. Solid crust of strange quark star 

The crust in our model is described by the BPS model of 
dense matter below the neutron drip (Baym et al. 1970). 
Neutrons are absorbed by SQM, and therefore the den- 
sity at the bottom of the crust, pb, cannot be higher than 
neutron drip density pnd — 4xl0^^g cm~^. A super- 
strong outward-directed electric field separates nuclei of 
the crust from the quark-plasma edge, which is necessary 
for the crust stability (Alcock et al. 1986). The gap be- 
tween nuclei and SQM surface should be sufficiently large 
to prevent the absorption of nuclei by SQM via quantum 
Coulomb barrier penetration. If at /Ond the gap is suffi- 
ciently large and the Coulomb barrier is sufficiently high, 
then pb = Pnd • Another possibility is that the gap and the 
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Coulomb barrier at pnd can be penetrated by nuclei, and 
then the stability conditions are reached at some lower 
value of pb- In such a case, both mass and thickness of the 
crust are lower than that corresponding to pb — Pnd (see, 
e.g., Huang & Lu 1997). In what follows, we will assume 
that pb = Pnd- The EOS of the crust is joined with that 
of the SQM core by requiring the continuity of pressure 
at the crust-SQM core interface. 



2.3. Calculation of models of stationary rotating 
strange stars with crust 

The general relativistic models of stationary rotating 
strange stars have been calculated by means of the 
code developed by Gourgoulhon et al. (1999), which re- 
lies on the multi-domain spectral method introduced by 
Bonazzola et al. (1998). With respect to the version 
of the code used in Gourgoulhon et al. (1999), Zdunik 
et al. (2000), Gondek-Rosihska et al. (2000), Zdunik & 
Gourgoulhon (2001), we introduce a second domain in- 
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side the star to describe its crust, resulting in a total of 
four domains. Using the same notation as in Sect. 3.2 of 
Gourgoulhon et al. (1999), these four domains are: Dl, 
which covers the strange quark interior; D2, which con- 
tains the crust; D3, which covers the strong field region 
outside the star (and where the ISCO is located) and 
D4, which extends to infinity, thanks to the compactifi- 
cation transformation u = 1/r. Following the technique 
explained in Bonazzola et al. (1998), the boundary be- 
tween domains Dl and D2 is forced to coincide with the 
transition between the strange quark matter and the crust, 
so that the huge density discontinuity between the two re- 
gions (cf. Fig. ^) is located at the boundary between the 
computational domains Dl and D2 and therefore does not 
give rise to spurious oscillations (Gibbs phenomenon) . The 
number of coefficients used in the spectral expansions (or 
equivalently number of grid points) is typically 33 in r 
and 17 in in each of the four domains. The global nu- 
merical error is evaluated by means of the virial identi- 
ties GRV2 (Bonazzola & Gourgoulhon 1994) and GRV3 
(Gourgoulhon & Bonazzola 1994), this latter being a rel- 
ativistic generalization of the classical virial theorem. For 
instance, for the configurations depicted in Figs. ^ and ^, 
the GRV2 (resp. GRV3) relative error is 6 ■ 10"^ (resp. 
4 ■ 10"''). The radius and the frequency of the marginally 
stable orbit, i?ms, and its frequency t'ms is determined in a 
standard way (see, e.g.. Cook et al. 1994, for the equations 
to be solved). 

Note that in all the results presented in this article, the 
equatorial radius of the star i?eq is defined as the circum- 
ferential radius (i.e. the length of the equator, given by 
the spacetime metric, divided by 2tt). In the nonrotating 
limit, this equatorial radius reduces to the surface value of 
the Schwarzschild coordinate R employed in the Tolman- 
Oppcnheimer-Volkoff equations. Note also that all non- 
equatorial coordinates (e.g. polar radius. Figs. § and ^ 
are obtained by rescaling of the radial quasi-isotropic co- 
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bare SS, i'„, = 



SS with crust, i'^^i=0 

SS with crust, i/„, = 641 Hz 




Fig. 1. Mass versus (circumferential) equatorial radius for 
nonrotating bare strange stars (SS, solid line), nonrotating 
SS with crust (dotted), and SS (bare and with crust) rotat- 
ing with the maximum observed frequency (641 Hz). The 
short-dashed line corresponds to a rotating bare SS and 
the long-dashed one to a rotating SS with crust (the open 
circle marks the point at which 641 Hz is the Keplerian 
frequency) . 

ordinate r (see the line element (6) in Gourgoulhon et al. 
1999) by the factor i?cq/?'cq- 

Let us consider a strange star, rotating at a frequency 
frot, with equatorial radius Req. If Rms > Req, then sta- 
ble orbits exist for r > i?ms; the ISCO then has the 
radius Rms and the frequency I'ms, and there is a gap 
of width i?i„s — -Req between the ISCO and the strange 
star surface. However, if R^s < Req, then i?isco = ^eq, 
'^isco — i^oih{Rcq)] an accretion disk then extends down 
to the strange star surface (or, more precisely, joins the 
stellar surface via a boundary layer). It should be noted 
that the recent analysis of the observational data strongly 
supports the existence of the marginally stable orbit in 
the case of 4U 1820-30 (Kaaret et al. 1999). 



3. Results 

3.1. Global stellar parameters 

In the present section we describe the global parameters of 
the rotating strange stars with a crust. We compare these 
results with those obtained for bare strange stars and with 
results of Glendenning & Weber (1992). 

In Fig. we present the mass versus radius relation 
for nonrotating strange stars and for stars rotating with 
a period equal to the smallest measured period of pulsars 
(1.56 ms, which corresponds to a rotation frequency of 
641 Hz). The last point on the strange star (SS) with crust 
curve corresponds to the Keplerian frequency of this star, 
equal to 641 Hz. Strange stars with smaller mass have 
thicker crusts and < 641 Hz. 

In the next figures we present global parameters of ro- 
tating strange stars (with and without crust) along some 
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Fig. 2. Rotation frequency versus angular momentum of 
SS along constant baryon number sequences. Each se- 
quence is labeled by the gravitational mass M of its non- 
rotating configuration. The last points of these curves cor- 
respond to the Keplerian configurations for bare SS: filled 
dots: Glendenning- Weber approximation of the rotating 
crust; open dots: exact calculations of the Keplerian fre- 
quency. 



evolutionary sequences with fixed total baryon numbers, 
equal to those of nonrotating configurations of given grav- 
itational masses (M = 1.2,1.4,1.6,1.75 Mq). Of course, 
due to the rotation, the actual gravitational mass of the 
star increases with increasing rotation frequency. It is the 
baryon mass Mb which is fixed for each sequence ("evo- 
lutionary sequence"). 

In these figures, the Keplerian configuration obtained 
in our exact calculations is marked by an open circle. For 
comparison, we also marked by filled dots the results of 
(Zdunik et al. 2000), obtained using the Glendenning & 
Weber (1992) model for the crust on a rotating strange 
star. The last point of each curve corresponds to the 
Keplerian configuration for the bare strange star. 

In Fig. H we present t'rot as a function of the angular 
momentum of the star J. We see that for a bare strange 
star the Keplerian configuration is not achieved via the 
increase of i/rot ■ On the contrary, close to this configuration 
there exists a region for which v^ot is a slightly decreasing 
function of J : this corresponds to spin up by the angular 
momentum loss. Therefore, v^f^ > t'rot(>/max) = t'K- 

This effect turns out to be suppressed by the existence 
of the crust. Although the crust does not change in prac- 
tice the curve Vrot(.J), the radius of a rotating star with 
a crust is significantly larger, setting the maximum value 
of J to be smaller by ~ 20% than for a bare strange star. 
However, due to the flatness of the function Vrot{J) for 
high J, the Keplerian frequency is nearly the same as in 
the case of bare SS. 

The reason why v^ot{J) does not increase close to Jmax 
is explained in Fig. ||. The increase of J results in flat- 
tening of the star (i?cq increases, so J increases without 




Fig. 3. The edge of the SQM core and the crust surface 
for a strange quark star with crust (baryon mass Mb = 
1.63 M©), rotating at v^ot = 121 Hz . The Keplerian fre- 
quency i^K — 1217 Hz (see Sect. 2.3 for the definition of 
the coordinates used in this plot). 

change in i^iot)- In other words, the Keplerian configu- 
ration is achieved due to the increase in the equatorial 
radius related to the deformation of the star. We see that 
the strange star with a crust is more oblate than bare SS 
for the same J. This is due to the fact that the crust is the 
most deformed part of the star. This effect is visualized in 
Figs. H and ^, in which we show the effect of rotation at 
Vrot = 1210 kHz on the structure of a strange quark star of 
Mb = 1.63 M0. As apparent in Fig. ^, the crust thickness 
at the equator is much greater than at the pole (about 
five times). Correspondingly, the density gradient in the 
polar ring of the crust is dramatically less steep than at 
the pole (see Fig. ^). Generally, the increase of i?oq/^poio 
with the total stellar angular momentum J depends on 
the baryon mass of the star: the lower the value of Mb, 
the steeper the increase of i?eq/^poio with increasing J. 
This is visualized in Fig. |l^. 

The global parameter which is relevant from the point 
of view of the stability of rapidly rotating stars with re- 
spect to spontaneous triaxial deformations is the ratio 
of kinetic to gravitational energy T/W (Shapiro & Zane 
1998). In Fig. ^ we have plotted T/W versus rotation fre- 



quency, for four evolutionary sequences of rotating strange 
stars with some fixed baryon number. The differences in 
the maximum value of T/W between the cases of bare 
strange stars and strange stars with crust are quite signif- 
icant. For strange stars with crust the maximum value 
T/W « 0.18 0.19, which is more than 20% smaller 
than in the case of bare strange stars, with maximum 
T/W « 0.25. 
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Fig. 4. Logarith m of density versus the radial coordinate 
(defined in Sect. 2.3) along the polar (dotted line) and 
equatorial (solid line) directions. The strange star model 
is the same as in Fig. ^. The shaded areas correspond to 
the crust. 
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M=1.6 Mg 
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M-1.2 Mg 



12 3 
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Fig. 5. Oblateness of the star presented as the ratio of the 
equatorial to polar radius, versus stellar angular momen- 
tum J, along some constant baryon number sequences. 
Each sequence is labeled by the gravitational mass M of 
its non-rotating configuration. Thin lines correspond to 
bare strange stars. 



3.2. Innermost stable circular orbit 

Previous calculations suggested that the gap between the 
ISCO and the stellar surface exists for Keplerian config- 
urations within a rather broad range of stellar masses 
(Stergioulas et al. 1999, Zdunik et al. 2000). Our calcu- 
lations, presented in Fig. indicate that for strange stars 
with crusts the gap exists only very close to the maxi- 
mum baryon mass. For lower masses (even 1.6 Mq for the 
SQMl model for which M^^^ = 1.8 Mq) there is no gap 
close to the Keplerian configuration . Thus if the observa- 
tional data support the existence of the gap, the theoreti- 
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Fig. 6. T/W ratio as a function of v^ot along some con- 
stant baryon number sequences. Notations as in Fig ^. 



cal lower limit on visco for rapidly rotating strange stars 
will be higher (asterisks in Fig. 

For bare strange stars the gap exists up to the 
Keplerian frequency. Even in the case of smaller mass stars 
for which < ^cq for intermediate frequencies (e.g. 
M — 1.2 Mq) the gap reappears at high rotation rates. In 
the extreme case of low mass bare strange stars (for which 
relativistic effects are negligible) the gap appears close to 
Keplerian frequency due to the large oblateness of the star 
(Zdunik & Gourgoulhon 2001). As a result, the frequency 
of the ISCO at the Keplerian limit is much smaller than 
for the nonrotating star (especially for low-mass stars). 

Of course, the smaller the mass, the larger the differ- 
ence between bare strange stars and strange stars with a 
crust, the size of the crust being a decreasing function of 
M . For low mass strange stars with crusts {M < 1.1 Mq) 
the gap between the surface and the marginally stable 
orbit vanishes independent of the rotation rate. In the in- 
termediate mass range (~ 1.1 1.3 M©) the gap exists 
for a nonrotating star, then disappears for larger rotation 
rates and reappears at high frequencies. In this case the 
difference between a bare strange star and a star with a 
crust is clearly visible (Fig. Due to the existence of 
the crust, the region of t'rot in which the ISCO is deter- 
mined by the stellar surface (i.e. i?ms < -Roq) is larger and 
thus the frequency t'isco is significantly smaller because 
the orbit of ISCO corresponds then to the equatorial ra- 
dius and not to the mass of the star (as in the case of 
the marginally stable orbit). As a result, the maximum 
possible value of the frequency of the particle in a stable 
circular orbit around the strange star with crust (which is 
determined by the point Rms = Rcq) is smaller by about 
100 Hz than for bare SS. 

For relatively massive strange stars with a crust (^ 
1.3 1.6 Mq), the gap exists for all rotation rates except 
those very close to i/k- The condition R^s = Rcq defines 
then the minimum ISCO frequency consistent with the 
constraint of the existence of the marginally stable orbit 
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Fig. 7. Orbital frequency of the ISCO versus the rota- 
tion frequency of strange stars, for the SQMl EOS of 
strange matter. Each curve corresponds to a fixed baryon 
mass, equal to that of a static strange star of gravita- 
tional mass indicated by a label. Along each curve, the 
angular momentum increases from J = (static configu- 
ration) to Jmax (Keplerian limit). Open circles correspond 
to Keplerian configurations of strange stars with crust 
(filled circles - approximation of the crust based on the 
results of Glendenning and Weber (1992)). Segments be- 
low the filled circles can be reached only by bare strange 
stars. Asterisks are defined by the condition i?nis = ^oq^ 
below these asterisks in the high frequency region and be- 
tween asterisks for 1.2 Mq case at intermediate rotation 
rate (thick part of the curve) i?,„s < R^q and the ISCO 
is determined by the stellar radius i?cq- Dash-dotted line 
corresponds to the bare strange star of 1.2 Mq. The dif- 
ference with respect to the star with crust in the region 
~ 100 — 1100 kHz is due to the fact that for bare strange 
stars the marginally stable orbit disappears at higher j^rot 
and reappears at lower Vy-ot leading to the visco larger 
than for SS with crust by ~ 100 Hz. 



1.3 kHz for 1.2 < M/Mq < 1.6 in the case of the 
SQMl EOS, asterisks in Fig.|). 

Only for stars with a mass close to the maximum 
for nonrotating configurations 1.6 1.8 Mq) is the 
marginally stable orbit located above the stellar surface 
at Keplerian rotation rate. 

As the Keplerian rotation corresponds to the minimum 
of t'lscOi the existence of the crust increases these values 
by 100-^200 Hz compared to bare strange stars. As one can 
expect, ;^sco is then equal to since ISCO is determined 
by the equatorial radius of the star rotating with Keplerian 
frequency. 

However it should be stressed that the lower limit on 
t^isco consistent with the assumption that a marginally 
stable orbit exists (i.e. R^s > Rcq) is significantly higher 
and is equal to ~ 1.27 kHz for SQMl model and 1.1 kHz 
for fine tuning of strange matter parameters, leading to 
the SQM2 EOS (Zdunik et al. 2000). For this last EOS 



^'isco < 1 kHz can be reached by bare SS for very rapid 
rotation. In this case the condition R^s > Rcq eliminates 
rapidly rotating models with ^'isco < 1-1 kHz, and the 
condition fisco =1.07 kHz with a non-zero gap can be 
fulfilled only by slowly rotating strange stars very close 
to maximum mass. Contrary to bare strange stars, when 
even for low mass stars rotating at the Keplerian limit, 
the frequency of the orbiting particle is limited by the 
marginally stable orbit and is approximately constant for 
a wide range of stellar masses (Zdunik & Gourgoulhon 
pool ), for strange stars with a crust the marginally stable 



orbit does not exist at the Keplerian limit and the in- 
nermost stable orbit coincides with the radius of the star 
(except for the mass close to Mmax), which rapidly in- 
creases as the mass decreases, leading to the decrease of 
the frequency on ISCO. 

3.3. Parameters of the crust 

Two parameters of the crust, relevant for our consider- 
ations, are its baryon mass A/b, crust £ind its equatorial 
thickness ^cqi defined in terms of the circumferential ra- 



dius introduced in Sect. 2.3. In the static case these quan- 



tities will be denoted as ^^^^^ and to, respectively. In 
the static case we will additionally consider crust contribu- 
tion to the stellar gravitational mass M^j.^^^. This quantity 
does not have a clear invariant meaning, reflecting the fact 
that gravitational energy cannot be localized in general 
relativity. Accordingly, M^^^^^^. should be considered only 
as the difference between the total gravitational mass M 
and the value at the interior/crust boundary of the met- 
ric coefficient m(R) ("enclosed mass") which enters the 
Oppenheimer-Volkoff equations. 

Let us start with the simplest case of a static strange 
star. For M > Mq, the ratio to/R is small, and this en- 
ables us to derive approximate expressions resulting from 
the Oppenheimer-Volkoff equations of hydrostatic equilib- 
rium 



^ 4nR^P, 

crust (.jyj 



2GM\ 
~R^) ' 



7 i?"c^ P, 



(1) 



7 - 1 GM phc^ ' 

where pb and Pb are the density and pressure at the bot- 
tom of the crust, and 7 is the (average) adiabatic index of 
the outer crust. 

For the BPS model (Baym et al. 1970) and pb equal to 
the neutron drip density, we have Pb = 7.8 10^^ dyn cm~^, 
Ph/Phc^ = 0.002, 7 = 1.28. This yields 



to = 0.65- 



^ ^6 



1 - 0.295— 1 km 

Re 



AC = 3.710-^^ (l- 0.295|^ 



Mr. 



(2) 



where Rq is the stellar radius in units 10^ cm and M, — 
M/Mq. The above formulae reproduce the exact results 
with an error lower than 6% in the case of to and 4% 
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Fig. 10. Same as Fig. H but for the thickness of the crust 



ness and mass of the crust of a nonrotating strange (equatorial and polar). Note that the two plots are shown 
star. Solid line: exact solution of the Oppenheimer-Volkoff using different scales and that for a nonrotating star equa- 



equation; dotted line: approximation given by formula 




Fig. 9. Baryon mass of the crust as a function of rotational 
frequency, along sequences of constant baryon number. 
The gravitational masses M of the nonrotating configu- 
ration of each sequence are 1.75, 1.6, 1.4, .1.2, IMq (from 
bottom to top) . Note that less massive stars have a larger 
crust (see also Fig m. 



for Mcrust for strange stars with mass between 1 Mq and 
2 Mq (maximum mass). A comparison between approxi- 
mate and exact values is displayed in Fig. ||. 

Rotation will increase both the mass and the equa- 
torial thickness of the crust. Results of Glendenning and 
Weber (1992) for teq as a function of frot (at fixed value 
of Mb) could be very precisely reproduced by a simple 
formula ioq(fiot) = • [1 + 0.7(i'rot/^'K)^], where vk is 
the Keplerian (mass shedding) frequency of the strange 
star (Zdunik et al. 2000). A similar formula reproduces 
very well the rotational increase of the baryon mass of the 
crust. 

Exact results for the crust on the static and rotating 
strange stars are presented in Figs. || and 10, The mass of 



torial and polar thicknesses are equal. 

the crust which can be supported by the rotating strange 
star can even be 2.5 times larger than in the case of a 
nonrotating configuration. Another very important effect 
is a strong oblateness of the star, resulting in large i?eq- 
As a result, the crust at the equator is much thicker than 
at the pole (see Figs. || and ||). 

The mass and thickness of the crust can be approxi- 
mated by a quadratic function of rotation frequency up to 
about 500 Hz. Thus the results of Glendenning and Weber 
(1992) describe very well the parameters of the crust of 
the star rotating with intermediate frequencies (say below 
one half of the Keplerian one). For higher values of i^rot, 
both crust mass and equatorial thickness increase rapidly 
and we have to use a polynomial including terms of much 
higher (sixth to eighth) order in Vj^ot- The approximate 
formulae are found to be 

A/B,crust = MB.crust(l + 0-24 ^.^^,^3 + 0.16 1^^,^^) , (3) 

icq(i^rot) = io • (1 + 0.4 1^1,^^ + 0.3 z.f„t,3) , (4) 

where — ^'rot/lO'^ Hz. Equation (^) works quite well 

nearly up to the Keplerian frequency, but the equato- 
rial radius increases very rapidly as one approaches the 
Keplerian configuration. In view of this, Eq. (^ can be 
safely used below 1 kHz (i.e. about 80% of the Keplerian 
value), where the change of the equatorial thickness is 
~ 100%. The difference between our quadratic approx- 
imate formula and the fitting results of Glendenning & 
Weber (1992) is a few percent. 

Another interesting effect of rotation on the crust 
structure deserves comment. As Glendenning & Weber 
(1992) pointed out, the polar thickness also increases (al- 
beit weakly) due to the decrease in the gravitational force 



8 



J.L. Zdunik et al.: The crust of rotating strange stars 



(resulting from the oblateness of mass distribution) at the 
pole. This effect is confirmed by our calculations (see bot- 
tom panel of Fig. |l^ 

4. Deconfinement heating in spinning-down 
strange quark stars with crust 

The spin-down of a magnetized strange star implies a de- 
crease of its oblateness. This process is particularly pro- 
nounced in the strange star crust. This change of shape is 
opposed by the elastic stress build-up in the Coulomb lat- 
tice of nuclei. The crust is relatively "soft" with respect to 
shear deformations, because its shear modulus constitutes 
only a few percent of the compression modulus (Haensel 
1997, Haensel 2001). Therefore, in a zeroth approxima- 
tion, the shape of the shrinking crust of a spinning-down 
strange star can be modeled by the equilibrium shape of 
a fluid. The spin-down implies some compression of the 
bottom matter elements of the crust, most pronounced in 
the equatorial belt. Consequently, as soon as the density 
exceeds the critical value pb, the bottom crust elements 
are engulfed by the SQM core. We note that if pb = Pnd 
then the absorption of crust matter is triggered by the 
neutron drip out of nuclei, while for pb < Pnd the absorp- 
tion proceeds via quantum tunneling of nuclei through the 
Coul omb barrier, separating them from the SQM core (see 
Sect. U). 

A typical spin-down rate for a relatively young 1000 yr 
old pulsar (like Crab) is P ~ 10~^^ s/s, where P is the 
pulsar period. Using the relation J^roti^rot — —P/P^, and 
expression (^ for the instantaneous baryon mass of the 
crust, we calculate the rate at which the baryon mass of 
the crust decreases due to the strange star spin-down: 



B. crust 



-9.6 10 



IPI 



14 



10-5 M, 







P 



10 ms 



10 



-13 



s/s 



gs 



(5) 



After crossing the SQM surface, nucleons (those in nu- 
clei and dripped neutrons) dissolve into u and d quarks 
(via strong interactions), which will be followed by the 
strangeness-changing equilibration process u + d — > s + u 
via weak interaction. The whole process is expected to be 
strongly exothermic, with an energy release per unit ab- 
sorbed baryon number q^cc — 10 — 40 MeV, its specific 
value depending on the assumed SQM model (Farhi and 
Jaffe 1984). Using Eq. (^), we estimate the total heating 
rate at the crust-SQM interface as 



9 10' 



33 



P 



10 ms 



-3 



IPI 



-^^B, crust 9dc 



10-I3s/S 

-1 



10-5 MfT. 10 MeV 



ergs 



(6) 



The typical deconfinement heating, estimated above, may 
therefore exceed the neutrino losses (luminosity) from the 
SQM core of a 1000 years old (or older) strange star with 
a crust, - lO^i (P/10 km)^ (T/10^ K)^ erg/s, and 



therefore could significantly influence (namely - delay) the 
cooling of this compact object (Yuan & Zhang 1999). We 
mention that the deconfinement heating, described above, 
is much more efficient than the roto-chemical heating, re- 
sulting from deviations of the composition of the SQM 
core of a spinning-down strange star from beta equilib- 
rium and studied in (Cheng & Dai 1996). 

5. Summary and conclusion 

We calculated exact 2-D models of stationary rotating 
strange stars of masses 1.2 — 1.75 Mq, covered with a 
crust of normal matter, up to the Keplerian limit. Despite 
its low mass 10"^ Mq) and small thickness (a few 
hundred meters), the solid crust plays an important role 
for these hypothetical objects. We considered sequences of 
uniformly rotating strange stars of a fixed baryon mass, 
and increasing values of angular momentum J, and we 
found that due to the large deformation near the Keplerian 
limit, the decrease of the maximum frequency, compared 
to that for bare strange stars, is very small. Close to the 
Keplerian limit, the stellar moment of inertia increases 
rapidly with rotation frequency, and this results in a weak 
dependence of the maximum frequency of uniform rota- 
tion, t'max7 on the value of Jmax, which in the presence 
of the crust is lower by as much as twenty percent com- 
pared to that of bare strange stars. On the contrary, the 
effect of the presence of the crust on the properties of the 
innermost stable circular stable orbit is important, and 
consists of about fifteen percent increase in t'isco at the 
Keplerian limit, as well as the disappearance of the gap 
between the ISCO and stellar surface for a broad range 
of stellar masses and rotation frequencies. The solid crust 
also lowers by as much as twenty percent the maximum 
T/W ratio which can be reached by a rotating strange 
star, from 0.23 - 0.25 to 0.18 - 0.20. 

Both the mass of the crust and its equatorial thickness 
increase with the rotation frequency. We have shown that 
this increase is very well approximated by simple and pre- 
cise analytical formulae up to 80% of the Keplerian fre- 
quency. For a rotation frequency below one third of the 
maximum, the increase with respect to the static case is 
quadratic in the rotation frequency. However, close to the 
Keplerian limit, the dependence is much more rapid, which 
is consistent with the characteristic behavior of Jivrot), 
mentioned in the previous paragraph. 
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